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Abstract 

In this paper we will extend a p-adic algorithm of Boyd in order to study the 
size of the set: 



My) 



n 



: — = (modp) 



.=1 ' 



Suppose that p is one of the first 100 odd primes and y G {1, 2, . . . ,p — 1}, then our 
calculations prove that \Jp{y)\ < oo in 24240 out of 24578 possible cases. Among 
other results we show that |Ji3(9)| = 18763. The paper concludes by discussing 
some possible applications of our method to sums involving Fibonacci numbers. 

1 Introduction 

The goal of this work is to extend an algorithm from David Boyd's paper, 'A p-adic 
study of the partial sums of the harmonic series". In particular, Boyd developed an 
algorithm which enabled him to calculate complete solution sets of the congruence 

n ^ 

^_=0(modp), (1.1) 

for 497 of the first 500 primes (his calculations did not finish for p G {83, 127,397}, 
and he accidentally omitted p = 509 from his table). Boyd's computations strongly 
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supported the hypothesis he drew from probabihstic models, namely that Eq. (jl.ip 
should only have a finite number of solutions in n for every prime p [1]. 
In this paper we will consider the much larger class of sums defined by 

n j 

We will extend Boyd's computational method to find complete solutions sets of 
the congruence Gn{y) = O(modp) for every prime p < 550, and for every p-adic 
integer y with = 0, with 338 exceptions where our calculations did not finish. 
Our computations rely upon results derived in Theorems 12.31 and 12.41 Table 1 
lists the number of solutions of the congruence = (modp) for p < 50 and 

yE{l,2,...,p-l}. 

This topic is also closely related to several unsolved problems in classical number 
theory. Boyd mentioned the connection between Eq. (jl.ip . Bernoulli numbers, and 
regular primes. It turns out that the zeros of Gn{y) (modp) are also related to 
the Wieferich primes. The elementary congruence Gp_i(2) = (modp), implies 
that Gp_i(2) = (modp) holds if and only if 2^~^ = 1 (modp^) . Primes that satisfy 
2^^^ = 1 (modp'^) are called Weiferich primes, and the only known examples are 
p = 1093 and p = 3511 (see [2] or [1]). Despite the fact that just two Weiferich 
primes have been discovered, the proof that infinitely non-Wieferich primes exist 
depends upon the ABC conjecture. Thus it would be extremely interesting to find 
an unconditional proof that "|{n : G„(2) = O(modp)}| = 0" holds infinitely often, 
since this condition implies that p is not a Wieferich prime. Our computations 
(and heuristics) seem to suggest that approximately 36% of primes satisfy this last 
condition. In general we believe that \{n : Gn{y) = O(modp)| < oo whenever 
Vp{y) = 0, and our calculations have confirmed that this is true in at least 24240 
out of 24578 cases for p < 550. 

Finally, we will conclude the paper with a brief discussion of some further pos- 
sible extensions of our method. For example, we will show that Boyd's algorithm 
can be applied to study the congruence 

n p 

^^ = 0(modp), (1.3) 
i=i 

when the -Fj's are Fibonacci numbers. It seems likely that a large number of inter- 
esting congruences similar to Eq. ()1.3I) can also be studied using our method. 

2 Elementary properties of JN(y) cind Gn(y) 

Although the primary goal of this paper is to determine complete solution sets of 
the congruence G„(y) = (modp), we can easily consider the more general case of 
JAr(y) when N e Z. 
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Definition 2.1. Let JNiv) be defined by 

JN{y) = {n:Gniy)=OimodN)}. (2.1) 

Before performing any calculations we will use elementary number theory to 
narrow the scope of our investigation. First notice that Gn{y) = (mod A^) if and 
only if Gn{y) = 0{inodp^) for every prime power dividing N. Therefore it is 
obvious that 

n 

M'...p'n"(y) = r\'^pAy^' (2.2) 

i=l 

whenever the pj's are distinct primes. Likewise, it is clear that for any prime p we 
must have a sequence of inclusions: 

Jp{y) 5 Jp2{y) 5 Jp3(y) . . . 

In Section 12.21 we will show that easiest way to calculate Jps (y) is to first determine 
Jp{y), and then to check whether or not Gn (y) = O(modp'^) for every n G Jp{y)- 
Surprisingly, it is often more difficult to determine Jp^iy) than Jp{y). 

Now we will discuss which values of y € Q need to be considered. Elementary 
number theory shows that calculating Jp^iy) is easy if fp(y) > because almost 
all of the terms in Gn{y) vanish modulo p^. Likewise, it is clear that Jp!>{y) = if 
fp(y) < 0. Finally, when yp{y) = we can appeal to the following proposition: 

Proposition 2.2. Suppose that i^p{y) = 0. If y = y {modp^), then 

Gn{y) = Gn{y) {modp') , (2.3) 

and it follows that 

Jps{y) = Jps{y). (2.4) 
Proof. Notice that Eq. (j2.3p is equivalent to the congruence 

y- y^ -y^ =o(modp^), 

i 

which follows trivially from the claim that for any j > 1 

=0(modp^). 

3 

To prove this claim, suppose j = j'p^ where (j',p) = 1. Since fp(y) = 0, rearrange- 
ment shows that 

(2//y)^^^"-l = 0(modp^+'). 

Since y/y = 1 (modp**), there exists a p-adic integer /3 such that y/y = 1 + j3p^ , and 
therefore 

(1 + p'^Y''^' -1 = (modp^+^) . 
Induction on 7 verifies this last equality. ■ 
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It follows immediately from Proposition l2.2l that Jp{y) £ { Jp(l), Jp{2), . . . , Jp{p — 1)} 
whenever i^piy) = 0. Recall that Boyd was able to calculate | Jp(l)| for p < 550 be- 
cause the set Jp(l) possesses a tree structure. In particular, an integer n can only 
belong to Jp(l) if the integer part of n/p also belongs to Jp(l). If the harmonic 
series is defined by 

Hn:=Gn{l)=y2-^ (2-5) 
7 

then for A; G {0, 1, . . . ,]5 — 1} and n > 

Hpn+k ^ — + Hk (modp) . (2.6) 
p 

It follows from Eq. (j2.6p that Hpn+k = (modp) can only hold if Hn = (modp) 
holds as well, thus the set Jp(l) inherits its tree structure from Eq. (|2.6|) . Boyd 
calculated Jp(l) using the following algorithm: 

Boyd's Algorithm: 



First check if Hk = (modp) for every k £ {1, . . . ,p — 1}. For each k that 
satisfies = O(modp), check if any elements of {Hp^, Hp^+i, Hp^+p-i} 
also vanish modulo p. Iterate this argument whenever Hp^+j = O(modp) 
for some j S {0, 1, ... ,p — 1}. The algorithm terminates in a finite amount 
of time if | Jp(l)| < oo. 



In order to apply Boyd's algorithm to the problem of calculating Jp{y), we will 
first need to prove that the set Jp{y) has a tree structure for y G {2, 3, ... p — 1}. 
In particular the the next theorem proves that Gn{y) satisfies a recurrence relation 
which reduces to Eq. (|2.6|) whenever y = 1. 

Theorem 2.3. Suppose that yp{y) = and k G {0, 1, . . . ,p — 1}, then 

Gpn+k {y") = + y'^Gtiy) + {'"—^^ Gp.i{y) (modp) . (2.7) 

Proof. First observe that for k € {0, 1, ... ,p — 1} 

Gpn+k [yn = z^ — + z^i^ — — — + — - — • 

Reducing modulo p, this becomes 
Gn (y""' 



Gpn+k {yn ^ — + - — r ^p-i(y) + ^"^^^(2^) (^od?') 



G„(vP) GniyP ) -GniyP) 

+ + ( ) Gp-i{y) + y"Gfc(y) (modp) 



p p V y - 1 
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Since Eq. (j2.3p shows that Gn yy^ j —Gn (y^) = (modp^), it is easy to check that 
Gn (y^^l ~ Gn {y^)i /p = (modp), and the theorem fohows.B 



Although Eq. ()2.7p does not apply to Gn{y) directly, we can still use it to 
calculate all of the integers belonging to Jp{y). In particular, Eq. (|2.3|) shows that 
whenever i'p{y) = 0: 

Jp{yn = Jp{y)- (2.8) 

Since Eq. (j2.7p shows that Jp (y^) has a tree structure, we can calculate Jp (y^) by 
simply modifying Boyd's algorithm to replace Hn with G„ (y^). 

2.1 A p-adic expansion for Gpn (y^) — ^"^^^^ 

Although in principle Jp{y) can be calculated by combining Boyd's algorithm with 
equations (j2.3p and (j2.8p . such a naive approach will rapidly exhaust a computer's 
memory. Boyd encountered a similar problem since both the numerator and de- 
nominator of Hn grow exponentially as functions of n. In order to calculate Hn 
efficiently, he used a p-adic series for Hpn — Hn/p- Boyd showed that for s > 2 there 
exists a polynomial, X{n), of degree s — 1 such that 

Hpn- — ^ X{n) (modp^) . 
p 

Notice that the right-hand side of this last congruence is easy to calculate even 
for extraordinarily large values of n. Thus if Hn (modp'^"*'^) is known, the value of 
Hpn (modp'^) follows easily from rearranging the congruence to obtain 

Hpn =X{n) + ^ (mod/) . 
p 

This argument can be iterated s — 1 times. For example, it is clear that 

H„2n = X{pn) + ^ (mod/-i) 
p 

- Yi ^ , Hpn{mo<lp') , 

= X (pn) H I mod p j 

p 

^X(pn) + ^ + ^(mod/-i) 

Although every iteration of this argument causes the loss of one digit of p-adic 
precision, we can still compute the value of Hn (modp) for every integer n £ Jp(l) 
by simply starting with a large enough initial value of s (assuming of course that 
|Jp(l)| < do). 

Therefore it is obvious that we will need to derive ap-adic expansion for Gpn {y^) — 
Gn (y^) Ip- Recall that Boyd determined his p-adic expansion for Hpn— Hn/p by first 
calculating the value of the function for n G {1, . . . s}, and then by multiplying those 



5 



values times the inverse of an s-dimensional Vandermonde matrix. Unfortunately 
Boyd's approach usually fails in this the numerator of Gpn (y^) — Gn {y^) /p 

grows far too quickly to allow for direct calculation. As a result, we will avoid those 
computations altogether by finding explicit formulas for the p-adic series coefficients. 
The first step is to split the function into two pieces: 

(i,p)=i 

In the next theorem we prove two useful formulas for calculating Eq. (j2.ip . 



Theorem 2.4. Suppose thatp > 2, = 0, andy ^ 1, then we have the following 

expansions: 

pn pj s-l 

^ y-^Aois)-yP'^^AisW{modp'), (2.2) 

j = l 2 = 

(i,p)=i 

E y'^'y"' ^No{s) - yP'^ Yl ^^(")"' (™«'^^'') • (2.3) 
j=i i=0 

We can calculate Ai{s) and Ni{s) using 

.(.^E(-.>KT)(|:^)(e«l''^) 



(modp^) , 

(2.4) 



N,{s)^-p'^^' — 



(i + 1)! 



where z is a p-adic integer defined by 

1 ^ (_i)i+i 



^ - 2 ^ 



and the 's are Stirling numbers of the second kind J^. 



Proof. We will prove equations (|2.2p and (j2.4p first. Rearranging the sum 
shows that 

3=1 •' k=l k=l r=l ^ 

a,p)=i 
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Reducing modulo this becomes 

j=l k=l k=l r=l \ \ k J / 

(i,p)=i 

Employing a geometric series yields 

y I \ y 1 i\ / 1^m„rr^ I \ » 1 I \ „m_.v r 

(J,P)=1 



j=l \k=l / m=0 \fc=l / \r=l 



(2.7) 



Now we will assume that y 7^ 1, then by Eq. (j4.2p 

n— 1 m ^ \ 

r=l j=0 ^ ^ 

Substituting this result into Eq. (j2.7p yields 

pn • s-1 



5^ ^^^o(5)-y^"2]^.(^K(mod/), 

i=i i=0 
(i,p)=i 

where 

aw = (-i)v(e|^) + E(-i)>"(7) (e^Ill 

\fe=l / m=i ^ ^ \A:=1 / 



Eq. (12. 4p follows from combining this definition of Aj(s) with Eq. ()4.1I) . 

Now we will prove equations (j2.3p and (j2.5p . First observe that if a = — l) /p^, 

then 

y:^^^^ = E^((i+V)'- 

i=i ^"^ i=i 
Applying the binomial formula shows that 

j = l j = l k = l ^ ^ 

= E^Et«^-^'--<^'-'+^»- 

k=i ■ j=i 
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Reducing this last equation modulo is easy. The nested sum is an integer, and 

k 2k— 1 

since i^p{y) = elementary number theory shows that t'p(^^-^ri — ) > k for all k. 
Therefore we can truncate the right-hand sum after the first s — 1 terms to obtain 

We will simplify Eq. (|2.8|) by using properties of the Stirling numbers of the first 
kind. Recall that for A; > 1 the Stirling numbers of the first kind have the generating 
function 

k 



^(x - 1) . . . (x - A; + 1) = ^ 4 



m=l 



Substituting this definition into the nested sum in Eq. (|2.8p yields 

k k-i 



E f ^^^^ - 1) • • • - ^ + = E T ^ ''^'^"^ = ^ s'^"''J:ry 

j=l i=l •' m=l m=0 j=l 



and therefore Eq. (|2.8p becomes 

E -EE sl""'^ Y.rr I (modp.) . (2.9) 

i=l m=0 \k=m+l ' / \j=l 

Now we will use a second power series identity for Stirling numbers. It is well 
known [3j that if |x| < 1 



1 (^v(-i)^"^.r_ v^M^^ 

m\\^ k ~ ^ k\- 

\k=l / k=m 



ml \ ^ k ] ^ k\ 

\k=\ / k=m 

It follows that for some integer- valued polynomial Q{x) 

\fc=l / k=m ^ 

where dg is the least common multiple of all integers less than s that are relatively 
prime to p. Taking x = ap and assuming that m > 1 and p > 2, we have 

Q (ap) > f p j- > s + m > s. 



mlcPp J \ ml J p — 1 

It follows that 

j=m ■ \j=l •' 
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and dividing both sides by p and then simphfying yields 

-z'^(mod/), (2.10) 



7! m! 

where z is defined in Eq. ()2.6p . Substituting Eq. (j2.10p into Eq. (j2.9p yields 

j=i m=o V ^ 7- y^-^i y 

Finally if y / 1 we can substitute equations ()4.ip and (j4.2p to complete the proof. 



2.2 Summary of computations 

In summary, we calculated at least part of Jp (y^) for every prime p < 550, and 
for each integer y S {2, 3, ... ,p — 1}. We then used the relation Jp (y^) = Jp{y) 
to determine Jp{y). We also checked Boyd's calculations of Jp(l) with a version of 
our program. In particular, we used Theorem 12.41 to determine Jp((l +pY), and 
then we verified Boyd's results from the fact that Jp((l + pY) = Jp{^)- Table 1 
lists the values of \ Jp{y)\ for p < 50, and an extended list for p < 550 is available at 
[www .math.ubc . ca/'^atrogers/Papers/padic .html, 

Some interesting observations follow from our computations. Firstly, \Jp{y)\ is 
small for many values of y and p. For example, when we considered the possible 
values of |Jp(y)| for p < 50, we found that only 27 out of 313 possible cases have 
\Jp{y)\ > 50. The three largest sets for p < 50 are | J47(12)| = 40608, | J47(8)| = 
27024, and |Ji3(9)| = 18763. In particular, we have explicitly proven that the 
congruence 

^ — = (mod 13) , 
i=i 

has exactly 18763 solutions. The first solution occurs at n = 3, while the largest 
solution has 419 digits and approximately equals n ~ 2.385 x 10^^^. 

We have also calculated that \Jp{y)\ = in 104 out of 313 possible cases for 
p < 50. This seems to agree with a simple heuristic suggesting that the density of 
such Jp(t/)'s should approach 1/e .36. To see this fact, notice that \Jp{y)\ = 
if and only if Gn{y) ^ (modp) for every n € {1, . . . ,p — 1}. If we assume that 
the value of Gn{y) (modp) is randomly distributed whenever y / 1 (recall that 
|<^p(l)| > 3 for any odd prime p), then it is clear that \ Jp{y)\ = with probability 
(1 — 1/pY~^ when y 7^ 1, and probability zero when y = 1. Therefore the expected 
ratio of empty Jp(y)'s for p < n equals 

E {\Jp{y)\ =0:p<n,l<y<p-l) = = _ -, 
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and the expectation approaches 1/e as n ^ oo by standard analysis. 
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13 
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19 
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29 


31 


37 


41 


43 


47 


1 





3 


3 


13 


638 


3 


3 


19 


3 


18 


26 


15 


3 


27 


11 


2 







37 








1 





9 


3 


2 


1 





29 








3 






4 


4 


184 





4 


4 
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140 
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1 





1 


5 


3 





1 





10 





5 








5 








12 


4 


1 





1 





1 
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34 


6 
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1 


16 





2 


6 


1 
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4 


4 





1 








129 
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4 


5 


1 


6 





325 


7 





27024 


9 










26 


18763 


1 


25 





1 





1 





4 


6 


10 










1 


2 


6 


1 





1 


1225 


27 


2 





1 


11 












6 





154 


14 


3 


2 


1 








4 


12 












11 


45 


13 





3 








17 





40608 


13 

















1 


2 








4 





1 





14 














1 





1 


2 


1 


133 


3 


13 


349 


15 














10 


4 


86 


1 


3 





1 


2 


1 


16 














65 


61 
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24 
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17 
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18 
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59 
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19 
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38 








20 
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1 
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151 
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5 





21 
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29 
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23 




















28 


48 





85 
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24 























24 


233 
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20 


92 
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10 


68 
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28 
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28 
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14 
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18 
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34 
























24 
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22 
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37 
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14 
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38 
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40 
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21 
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Table 1 : Values of \Jp{y)\ for p < 50 
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7 


11 


13 


17 


19 


23 


29 


31 


37 


41 


43 


47 


43 

44 
45 
46 






























6 
1 
5 
2 



Table 1 (continued): Values of | Jp(y)| for p < 50 



We can examine cases where y is fixed and p varies in somewhat greater detail. 
Following Boyd we will use the notation 

Gm{y) = Jp{y) [^{p'^-\p'^-^ + 1, . . . - 1}, 

to denote a level in the tree Jp{y)- Notice that the integers contained in Gm{y) are 
precisely the elements of Jp{y) that we will discover during the m'th iteration of 
Boyd's algorithm. As usual, Mp{y) — 1 equals the number of levels in the tree Jp{y), 
notice that 

Mp(,y)-l 

My)= U ^rniy). 
m=l 

As an example we will consider the case that occurs when y = 2. The following 
table lists the nonzero values of | Jp(2)| for p < 550. 



p 


Mp(2) 


l^p(2)| 


Values of |Gm(2)| for 1 < m < Mp{2) 


5 


15 


37 


1,2,3,4,2,3,3,4,3,2,4,1,2,3 


13 


2 


1 


1 


19 


5 


9 


1,2,3,3 


23 


3 


3 


1,2 


29 


3 


2 


1,1 


31 


2 


1 


1 


41 


9 


29 


1,3,3,5,7,7,2,1 


53 


2 


2 


2 


59 


5 


7 


2,2,1,2 


73 


6 


11 


3,2,1,2,3 


83 


3 


2 


1,1 


89 


24 


56 


1, 1, 1, 2, 4, 4, 7, 6, 4, 3, 2, 2, 2, 2, 1, 1, 1, 3, 3, 3, 1, 1, 1 


103 


3 


3 


2,1 


113 


50 


394 


1, 3, 4, 7, 9, 10, 10, 9, 9, 7, 7, 11, 8, 7, 6, 7, 10, 12, 15, 11, 9, 12, 9, 7, 7, 
6, 11, 8, 8, 11, 14, 14, 14, 11, 11, 12, 9, 11, 6, 7, 3, 7, 6, 3, 5, 3, 2, 2, 3 


131 


18 


80 


2, 3, 4, 6, 2, 3, 3, 4, 3, 6, 3, 6, 12, 11, 7, 4, 1 


137 


6 


9 


1,2,2,3,1 


151 


3 


3 


2,1 


157 


11 


18 


2,1,1,3,1,2,2,3,2,1 


163 


2 


1 


1 



Table 2 : Primes p < 550 for which \Jp{2)\ > 0. 
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p 


Mp(2) 


l^p(2)| 


Values of |G,ri(2)| for 1 < m < Mp{2) 


167 


2 


1 


1 


173 


3 


3 


2,1 


179 


2 


1 


1 


181 


2 


1 


1 


193 


4 


3 


1,1,1 


197 


6 


7 


2,1,1,1,2 


199 


4 


3 


1,1,1 


211 


12 


41 


1,2,5,5,5,6,3,6,4,3,1 


239 


3 


7 


3,4 


241 


6 


9 


2,2,2,2,1 


257 


9 


17 


1,1,1,3,3,2,4,2 


269 


7 


7 


2,1,1,1,1,1 


271 


6 


7 


1,3,1,1,1 


293 


2 


1 


1 


307 


4 


6 


1,3,1,1 


311 


3 


3 


1,2 


313 


2 


2 


2 


317 


4 


7 


3,3,1 


331 


13 


55 


1,1,2,3,3,6,10,6,9,9,3,2 


337 


20 


47 


2,1,3,5,1,3,4,4,2,3,3,3,2,3,1,2,3,1,1 


349 


3 


4 


1,3 


367 


2 


1 


1 


373 


4 


5 


2,2,1 


379 


6 


19 


4,3,6,4,2 


383 


3 


2 


1,1 


389 


17 


51 


1,2,2,2,2,5,7,8,6,6,3,2,1,1,1,2 


397 


13 


33 


1,3,3,4,3,1,3,2,6,3,3,1 


401 


2 


1 


1 


419 


2 


1 


1 


431 


12 


76 


3,6,4,10,11,10,8,8,9,4,3 


439 


14 


26 


1,1,1,1,1,1,3,4,3,3,2,3,2 


449 


3 


2 


1,1 


457 


7 


7 


1,1,2,1,1,1 


461 


2 


1 


1 


463 


8 


12 


1,2,2,1,3,2,1 


479 


2 


2 


2 


487 


21 


52 


2,2,1,3,3,4,6,2,3,4,2,3,2,2,1,3,4,2,2,1 


499 


30 


272 


1,6, 6, 6, 6, 11, 9, 11, 10, 16, 15, 18, 14, 18, 16, 11, 
10,11,8,6,9,9,10,8,7,5,8,5,2 


509 


5 


4 


1,1,1,1 


523 


8 


16 


2,2,2,2,4,3,1 


547 


4 


4 


1,2,1 



Table 2 (continued) : Primes p < 550 for which | Jp(2)| > 0. 
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Finally, we will point out that it is usually easy to calculate Jps [yP) after first de- 
termining Jp{y). Since we will have already calculated the value of Gn (y^) ^modp'''^ 

for some s' ^ 1, we can simply check whether or not Gn (y^) = (modp^) for every 
n G Jp{y) and for any j < s'. In practice this check rarely requires new computa- 
tions, since generally we will have used a value of s' much larger than the order of 
vanishing of G„ (y^) modulo p. As an example we proved that 

J5(2) = {3, 17, 19, 86, 97, 99, 485, 488, 497, 499, 2486, 2496, 12431, 

12482, 12484, 62157, 62159, 62421, 310787, 310789, 312107, 
312109, 1553936, 1560537, 1560539, 7802685, 7802688, 
39013425, 39013428, 39013442, 39013444, 195067126, 
975335630, 975335633, 4876678152, 4876678154, 4876678166}, 

and with minimal extra computations we also determined that 

J25(7) ={3, 19, 499, 2486, 12431, 312107}, 
Ji25(32) =0. 

Notice that J^s [32) = when s > 3, since in those cases J5s(32) C Ji25(32) = 0. 

Unfortunately this procedure is unsuitable for calculating the majority of val- 
ues of Jps{x) when s > 1. Although we can easily calculate Jps (y^), the method 
only applies to Jps{x) when x = y^ (modp^) for some y. For example, if x G 
Z25\{1, 7, 18, 24}, then we have to settle for the weak conclusion that J25ix) C J^ix). 
While Theorem 12.41 probably only requires minor modifications to extend the com- 
putations, we will not address that problem here. 

Open Problem : Calculate Jpa{x) when s > 1 and x ^ y^ (modp*) for any y. 



3 Conclusion 

Although we primarily restricted our attention to computational problems in this 
paper, it would be desirable to construct probabilistic models to explain the behavior 
of \Jp{y)\. Boyd constructed such models to explain the behavior of |Jp(l)|, and it 
seems likely that his ideas will suffice to explain our results as well. Notice that 
after combining our calculations with Boyd's, we have proved that \Jp{y)\ < oo in 
98.6% of cases for primes p < 550. In fact, it would be desirable to find a general 
proof of the following conjecture: 

Conjecture 1 : We will conjecture that | Jp(y)| < oo whenever fp(y) = 0. 

Finally, we will conclude the paper with an example of a more complicated type 
of function that we can easily study with our method. Consider the function 
defined by 

" ^. 
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where the Fj^s are Fibonacci numbers. Recall that we can either calculate Fj re- 
cursively, or with Binet's formula: 

It is not difficult to prove that the solutions of /„ = (modp) are arranged in a 
tree. The crucial fact for proving this claim is that the Fibonacci numbers satisfy 
the congruence 

Fpsj = (I) Fps-i^ (modp^) , (3.1) 



implies that fn^ = (I) /« (modp), where 



for all s and j, with (^) denoting the Legendre symbol. This congruence easily 



f(i) _ V ^Pi 

n 2^ j ■ 



Using properties of the Fibonacci numbers we can prove that fn^ obeys the con- 
gruence 

P 



^ fc=i fc=i 

and therefore our claim about the distribution of zeros of /„ (modp) follows im- 
mediately. Based on cursory computations it also seems reasonable to make the 
following conjecture: 

Conjecture 2 : For all n we have f^n = (mod 5). Furthermore, if p 7^ 5: 

|{n : /„ = (modp)}| < 00. 

From these short computations, it seems obvious that our method will extend 
to functions involving integer sequences other than just the Fibonacci numbers. We 
will speculate that many functions of the form 

should obey p-adic recurrences analogous to equations (|2.4p or (j3.2p . provided that 
Tj satisfies a second-degree linear recurrence. We are not presently prepared to 
speculate on the behavior of such functions when the T^-'s satisfy higher order re- 
currences, as we failed to observe any interesting patterns modulo p when the Tj 's 
equal Tribonacci numbers. 
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4 Appendix : A simple but important sum 

Virtually all of the calculations in this paper depended upon our ability to efficiently 
calculate the simple sum 

n 

i=i 

for extremely large, hut finite^ values of n. It was therefore imperative to eliminate 
the n-dependency from the index of summation. While many obvious formulas exist 
for this sum, including 



we chose to avoid recursive identities, and to find a closed form instead. Perhaps 
the crucial observation was that when n = oo: 

oo ^ i 

Li_,(x) := ^/x^- = -5ro + Ei!6(^) \ . (4.1) 

As usual 5^0 is the Kronecker delta, and S^"''* denotes the Stirling numbers of the 
second kind [3j. Therefore, briefly assuming that |x| < 1, we obtain 

n oo oo 

i=i i=i i=i 

and expanding (n + j)** with the binomial formula yields 

^/x^ =Li_,(x)-x" ( nn™Li_(,„^)(x). (4.2) 

j = \ 7Tl=0 ^ ^ 

Since the left-hand side of this last identity is a polynomial, the principle of analytic 
continuation shows that the identity holds for all x. In practice we should only use 
Eq. ()4.2p if X 7^ 1, since in that case we can use Eq. (j4.ip to calculate Li_r(x). 
When X = 1 we can calculate the left-hand side of Eq. (j4.2p by simply reverting to 
Bernoulli's classical formula for power sums. 
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